Efficient algorithms for constraining orientation tensors in Galerkin methods for the Fokker-Planck equation

نویسنده

  • Christoph Lohmann
چکیده

This paper deals with the problem of positivity preservation in numerical algorithms for simulating fiber suspension flows. In contrast to fiber orientation models based on the Advani-Tucker evolution equations for even-order orientation tensors, the probability distribution function of fiber orientation is approximated using the Galerkin discretization of the Fokker-Planck equation with Fourier basis functions or spherical harmonics. This procedure leads to a natural generalization of orientation tensor models replacing ad hoc closure approximations by Galerkin equations for the fine-scale components. After introducing an operator splitting approach to solving the discretized Fokker-Planck equation, we present conditions and correction techniques that guarantee physically correct distribution functions. As the reader will see, the derivation of these conditions is independent of the space dimension and their applicability is not limited to the simulation of fiber suspensions.

برای دانلود متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید

ثبت نام

اگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید

منابع مشابه

Pseudo-spectral ‎M‎atrix and Normalized Grunwald Approximation for Numerical Solution of Time Fractional Fokker-Planck Equation

This paper presents a new numerical method to solve time fractional Fokker-Planck equation. The space dimension is discretized to the Gauss-Lobatto points, then we apply pseudo-spectral successive integration matrix for this dimension. This approach shows that with less number of points, we can approximate the solution with more accuracy. The numerical results of the examples are displayed.

متن کامل

Analysis and Implementation of Numerical Methods for Simulating Dilute Polymeric Fluids

Analysis and Implementation of Numerical Methods for Simulating Dilute Polymeric Fluids David Knezevic Doctor of Philosophy Balliol College Michaelmas Term 2008 In this thesis we develop, analyse and implement a number of numerical methods for simulating dilute polymeric fluids. We use a well-known model in which the polymeric fluid is represented by a suspension of dumbbells in a Newtonian sol...

متن کامل

A Fully Discrete Discontinuous Galerkin Method for Nonlinear Fractional Fokker-Planck Equation

The fractional Fokker-Planck equation is often used to characterize anomalous diffusion. In this paper, a fully discrete approximation for the nonlinear spatial fractional Fokker-Planck equation is given, where the discontinuous Galerkin finite element approach is utilized in time domain and the Galerkin finite element approach is utilized in spatial domain. The priori error estimate is derived...

متن کامل

Numerical Studies and Simulation of the Lower Hybrid Waves Current Drive by using Fokker – Planck Equation in NSST and HT-7 Tokamaks

Recent experiments on the spherical tokamak have discovered the conditions to create a powerful plasma and ensure easy shaping and amplification of stability, high bootstrap current and confinement energy. The spherical tours (ST) fusion energy development path is complementary to the tokamak burning plasma experiment such as NSTX and higher toroidal beta regimes and improves the design of a po...

متن کامل

ذخیره در منابع من


  با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید

برای دانلود متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید

ثبت نام

اگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید

عنوان ژورنال:
  • Computers & Mathematics with Applications

دوره 71  شماره 

صفحات  -

تاریخ انتشار 2016